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Abstract 

We propose an explicit expression for vacuum expectation values ( e ia<p ) of the expo- 
nential fields in the sine-Gordon model. Our expression agrees both with semi-classical re- 
sults in the sine-Gordon theory and with perturbative calculations in the Massive Thirring 
model. We use this expression to make new predictions about the large-distance asymp- 
totic form of the two-point correlation function in the XXZ spin chain. 
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Many 2D Quantum Field Theories (QFT) can be realized as Conformal Field Theories 
(CFT) perturbed by some relevant operators In these cases the correlation function of 
the local fields O a (x) of the perturbed theory 

(O ai (x 1 )...O aN (x N )) (1) 

are defined in a formal way through the Conformal Perturbation Theory (CPT) as 

Z^y 1 (O ai (x 1 )...O aN (x N ) e-»f d2x * (x) ) C FT, (2) 

where 

Z^) = {e-^ d2x ^ x) ) CFT , 

$ is the local field taken as the perturbation, \x is the CPT expansion parameter and 
{•••)cft denote the expectation values in the original CFT. Of course, the expression (0) 
for the correlation functions ([D can not be taken literally. As is well known, the power 
series expansion in \i generated by (H) is always plagued by infrared divergences, reflecting 
the fact that the correlation functions ([!]) contain non-analytic in \x terms. However these 
non-analytic terms are rather the property of the vacuum | ) in (...) = (0|...|0)/(0|0) 
then of the local operators O a . The local fields are believed to satisfy the OPE algebra 

O a (x) O b (y) ~ C ab( x ~ V) °c(y) ■ (3) 

c 

By the very virtue of these relations the c-number coefficient functions C are not very sen- 
sitive to the long-range properties of the theory (in particular, the CPT for the coefficients 
C do not contain any infrared divergences) and therefore one expects that they can be 
expanded into a power series in \x with finite radius of convergence. On the other hand, 
by successive application of (|3|) any correlation function (p can be reduced down to the 
one-point functions (O a (x) ). Simple dimensional analysis shows that these quantities are 
non-analytic in \x 

(O a (x)) = G a ^ , (4) 

where A a are conformal dimensions of the fields O a and G a are constants. In writing (]4|) 
we have assumed that the perturbing operator $ consists of a single scaling field of the 
conformal dimension A. In this sense it is the one-point functions (|]) that contain all the 
non-perturbative information about the QFT which can not be extracted directly from 
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CPT 0]. Therefore the problem of calculation of the one-point functions has fundamental 
significance 0. 

During the last two decades much progress has been made in understanding 2D In- 
tegrable QFT. However the closed expressions for the one-point functions are obtained in 
very few cases. Most significant result obtained in this direction so far was exact computa- 
tion in many cases the expectation value of the perturbing operator $ in (0) with the help 
of Thermodynamic Bethe Ansatz and similar techniques 0. In this paper we propose an 
explicit expression for the one-point functions of all exponential fields O a (x) = e iaip (x) in 
the sine-Gordon QFT. 

The sine-Gordon model is defined by the Euclidean action 

A SG = j d 2 x | y^(<5M 2 - 2// cos(/3^) J , (5) 

where (3 and \x are parameters @. The renormalization required with this definition is 
relatively simple. The constant (3 does not renormalize while \i is renormalized multiplica- 
tively. In order to give the constant \x a precise meaning one has to fix the normalization 
of the field cos(/3(p). A convenient normalization which we will accept here corresponds 
to the short distance limit of the two-point function 

1 2 

( cos(f3(p)(x) cos (/?</?) (y) ) — > - \x — y\~^ as \x — y\ — > . (6) 

This normalization is the most natural one if one understands (|5|) as a perturbed CFT, 
with the field cos(Pp) taken as the perturbing operator. Under this normalization the 
field cos(/3<£>) has the dimension [length ] and correspondingly the dimension of the 
constant \x is [length] 213 

The theory (|5|) has a discrete symmetry p — > (p + 2nn/[3 with any integer n. In this 
paper we study (j5|) in the domain < /3 2 < 1, where the above symmetry is spontaneously 
broken (and the QFT (|5|) is massive), so that the theory has infinitely many ground states 
| n ) characterized by the associated expectation values of the field (p, ((p) n = 2im//3, 
where ( . . . ) n = ( n | . . . | n )/( n | n ). In what follows we concentrate our attention on 
one of these ground states, say | Oq ), and use the notation (...) = ( . . .)o- 



1 Although our discussion explicitly concerns 2D QFT, the general statements above are ex- 
pected to hold in multidimensional QFT as well. 

2 Note that our notation j3 for the sine-Gordon coupling constant differs by a factor \Z8tt from 
conventional one 0. 
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The sine- Gordon model admits an equivalent description as the Massive Thirring 
model j| 

A M t = J d 2 x | i^d^-M^- 9 -^!^) 2 | , (7) 

where is the Dirac field, the four-fermion coupling constant g relates to (3 in (|5|) as 

9 1 



7T 2/3 2 



1, (8) 



and M. is the mass parameter. In particular, the fermion current of ([?]) is related to the 
field (p in (|5|) as 

j» = ^ = -JL f*d x <p. (9) 

The QFT (|5|) is integrable and its on-shell solution, i.e. spectrum of particles and 
S- matrix, is well known 0. The theory contains soliton S, antisoliton S (these particles 
coincide with the "fundamental fermions" of the Lagrangian ([?])) and a number of soliton- 
antisoliton bound states ( "breathers" ) B n , n = 1,2,...,< 1 /£ ; here and below the notation 

t=T^ (10) 

is used. The lightest of these bound states B\ coincides with the particle associated with 
the field ip in perturbative treatment of the QFT (|5|). Its mass m is given by 

m = 2M sin (tt£/2) , (11) 

where M is the soliton mass. Exact relation between the particle masses and the parameter 
H in ([D defined with respect to the normalization @ was found in 



r(/3 2 ) 

ii - 



M 



2 r(|) 



(12) 



This allows one to derived the vacuum expectation value of the field exp(ip(p) using an 
obvious relation 

<e^) = -^/(/x), (13) 
where f((i) is the specific free energy f(fi) = — y logZ, which is also known exactly 0, J7[ 

M 2 

/(//) = -— tan «/2) . (14) 
3 



Combining (|TI|)-([r3|) one finds 

(i + 0^r(i-/3 2 ) 



r(| + |)r(i-|; 



2/3 2 -2 



in 



2(3 2 



(15) 



16 sin(7r0 T(f3 2 ) \ 4^ 

where we have assumed that the field e 1 ^^ is normalized according to (|6|). 

In this paper we study the expectation values ([|) of the exponential fields O a (x) = 
e iaif {x) in the sine-Gordon model. Like in (|5|), we assume that these fields are normalized 
in accordance with the following short distance limiting form of their two-point functions 



e^(x)e-^(y)) 



\x - y\ 



-4a 2 



as 



\x - y\ 



. 



(16) 



Under this normalization the field e %a<p has the dimension [ length ] 
For two special values of /3 2 the one-point function 



-2a 1 



(17) 



admits direct calculation. 

a) The semi-classical limit f3 2 — > 0. In Appendix A we study the expectation values 
( [FTP with a = a/ (3 for (3 2 — > and fixed a, and obtain the result 

sinh 2 (2at) 



2 l - 

m \ 0* 



dt 



t, 



t sinh(2t) 



-2a 2 e~ 2t 



(18) 



where D(a) can be calculated as the functional determinant (A. 11). 

b) p 2 = 1/2. In this case the Thirring coupling constant (|8|) vanishes and the sine- 
Gordon model is equivalent to the free-fermion theory. In Appendix B we use this simpli- 
fication to obtain for (3 2 — 1/2 and 3?eo I < l/y/2 



( e ia<p ) 



2a z 



YJ exp 



00 dt_ 
t 



smhVgog 2a 2 c -2t 
sinh 2 (t) 



(19) 



where M is the mass of the free fermion field. 

Evident similarity between (|18D and (jrj) suggests the following expression for the 
expectation value (|T7j) for generic (3 2 < 1 and hftea < 1/(2/3) 



Qa 



m 



exp 



£(i+j)£(Lii) 

' dt " 



2a z 



X 



sinh 2 (2a/%) 



o t [ 2sinh(/3 2 t) sinh(t) cosh ((1 - (3 2 )t) 



2a 2 e 



-2t 



(20) 



The pre-exponential factor in (|20| ) is chosen to agree with (15). The formula (ECT) is our 
conjecture for (p~7|) . 

The expression (^Up can be expanded in power series in /3 2 or in g (^) and the co- 
efficients of these expansions can be compared with the perturbative calculations for the 
actions (|5j) and ([F]), respectively. To perform this check it is more convenient to consider 
the expectation values ( ip 2n ) which can be readily obtained by expanding fllTD in power 
series in a 2 . Let us define "fully connected" one-point functions am as 



fa *> = i + £ 



so that 



a 2 = ((^ 2 ); a 4 = (^ 4 )-3(^ 2 ) 2 
Expansion of (^Dp gives 

T(i + |)r(i-|) 



(21) 



(22) 



2 + 2 



m. 

a 2 = -4 log <^ 

00 eft 

o T I sinh(/3 2 t) sinh(t) cosh ((1 - (3 2 )t) 

n A 2n-1 n2n ' 1 UL 



(3 2 t 



— e 



-2t 



(23) 



t7 2n = (-l) n 4^- 1 J 9 



o sinh(/3 2 t) sinh(t) cosh ((1 - /? 2 )t) 



n > 1 



We show in Appendix C that a 2 and a 4 in ( |23|) agree with the perturbation theory for @ 
up to /3 4 , and that a 2 agrees with the perturbation theory for ([7j) up to g. 

Clearly, more checks of (|20| ) are desirable. We note in this connection that the ex- 
pectation value ( |i~7D controls both short and long distance asymptotics of the two-point 
correlation function 

G aa >(\x-y\) = (e ia v(x)e ia 'v(y)) (24) 

with |a + a'| < (3/2. Indeed, if this inequality is satisfied the short distance limit of ( f24"| ) is 
dominated by OPE 



e ia v(x)e ia ' v (y) -> \x - y\ 4aa ' e i{a+a ' )v (y) 



as 



I a: — y\ — > . 



(25) 



Therefore 



<3aa>{r) 



Q , , r 4aa ' 

£a Ga> 



as 
as 



r 
r 





• 00 



(26) 



It is probably possible to check this relation by calculating numerically the correlation 
function (|24f ) with the use of exact form-factors |§. It is also worth noting that the 
expression ( p0|) is expected to hold for the one-point functions of the sinh-Gordon model 
as well, if one makes the substitution /3 2 — ► — /3 2 in (|20|). The form- factors for the sinh- 
Gordon model can be found in ||, 0. 

The one-point functions (|I7p of the sine-Gordon model can be used to derive the one- 
point functions of primary fields in c < 1 "minimal" CFT perturbed with the operator 

$1,3- As is known, the perturbed "minimal model" Ai p / P ' (with c = 1 , ) can be 

obtained by "quantum group restrictions" from the sine-Gordon QFT (|j) with £ 



P'-p 



10|j , |l 1|] . Using this relation one obtains from ( P0|) 



<*i,*> 



M 



Q({£ + l)l-£k) , (27) 



2 r(|) 

where the function Q(rj) for 1 9fte r\ | < £ is given by the integral 

dt ( cosh(2t) sinh(t(?7 - 1)) sinh(t(r} + 1)) (r? 2 - 1) _ 4t 



Qfa) ex P(_/ o t y 2cosh(t)sinh(tf)sinh(t(l+0) 2 £(£ + 1) 6 

and it is defined through analytic continuation outside this domain. In (|27j ) (J?^ stands 
for the primary field of the dimension 

_ ((g + i)f-gfe) 2 -i 

Ai '' i?«TT) (28) 

with canonical normalization 

($i,k(x) *iM) -> l^-2/r 4A '' fe as |x-2/|^0, (29) 

and M denotes physical mass in the perturbed theory. 

Finally, we remark that accepting the conjecture (|2(]) one can make an interesting 
prediction about long distance asymptotic of two-point correlation function in the XXZ 
spin chain 

F(f3 2 ) n' 02 as n -> oo , (30) 



( vac 




vac ) 


(vac 


vac 


) 



where | vac ) is the ground state of the XXZ Hamiltonian i 

l-(3 2 
2 sin(7r/? 2 



1/32 °° 

^xxz = . £^ Yl + + cos(7r/3 2 )(a s Vf +1 - 1)) , (31) 



3 As usual, the infinite XXZ spin chain is defined as N — > oo limit of a finite chain of size N 
with periodic boundary conditions. 
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a x s , a v s and a z s are Pauli matrices associated with the site s of the chain. While the power- 
like behavior in (130) is well known, the factor 



:i+o s 



r(f) 



X 



exp 



2^Fr(| + |) 

00 (ft / sinh(/3 2 t) 



/ t V s inh(i)cosh((l-/3 2 )t) 
in ( PP[) is the consequence of ( POD (see Appendix D). 



/3 2 e" 2t ) 



(32) 
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Appendix A. 

The expectation value flTTj ) can be calculated as the Euclidean functional integral 

Q a = Z- 1 ^) J [XV] e ia ^°> e~ AsG , (A.l) 

where Asg is the action (|5|). For a = a/ [3, (3 — > this integral is dominated by a saddle- 
point configuration <p(x) = <f>(x)/j3 which solves the classical equation 

dl <p(x) = m 2 sin <j>{x) - 8nia 5 {2) (x) , (A.2) 

and decays sufficiently fast at \x\ — > 00. Obviously, the solution depends only on the radial 
coordinate <fi = ^( r ) 5 r = \ x \- Note that is real only if a is imaginary, a = iuj. The 
classical action calculated on this solution diverges. Proper treatment of this singularity 
requires introducing small cutoff distance e and taking the limit 

log(e-^ (0) >~-^ lim (u 2 log(e 2 )-u<j ) (e)-2S(u 1 e)) , (A.3) 
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where 



S(u,e) 



rdr 
"16" 



(d r </>) 2 + 2m 2 (l-cos. 



and 4>{r) is the solution to the equation 



d^cj) + r 1 d r (f> — m 2 sin = 



(A.4) 



(A.5) 



which satisfies the asymptotic conditions 



4>(r) 



/mr\ ^ 



as r — > 



(A.6) 



0(r) 



■— sinh(7ro;) Xo(mr) as r — > oo . 



7T 



(A.7) 



Here C w is certain constant and Ko(t) is the Macdonald function. Let us introduce the 
function 

/ _ /mr\ \ 

(A.8) 



Taking derivative of ( |A.3| ) with respect to uj and using the equation (|A.5| ) one can derive 
the relation 

,., /„s 9 / _ /m\ f At \ 

(A.9) 



log( e-t"<» > ~ A ( - c 2 log (| ) + w / $ 5(r; 



As is shown in Appendix B the function ( A.8|) admits the following representation 



S(u) = uf 



00 dv ( cosh 2 7r(^ — u>) cosh 2 7r(z^ + u)) 

l°g o 

27T \cosh ttu cosh7r(z/ — 2u) cosh7r(z/ + 2u) 



MeV(--iv) , (A.10) 
2 



where *(t) = r'(t)/r(t). The exponential factor in (|T|) is derived from ([OP and (|OC| ) 
with the use of the integral representation for The pre-exponential factor in (|I8| ) 

can be obtained by evaluating the functional integral ( |A.1|) in the Gaussian approximation 
around the above classical solution d> 



2a/ 



Det' 



-<9 2 + ?Tl 2 COS ( 

—d 2 + m 2 



(A.11) 



where 7 = 0.577216... is the Euler constant. The first factor in ( |A.11|) appeared because 
of the mass renormalization in ( |A.9| ); its main role is to cancel ultraviolet divergences in 
the determinant. Note that validity of (^0|) requires the determinant ( |A.11|) to be 



D(iuj) 



exp 



1 f 00 dt sin 2 (2(x;t) 

2 Jo t cosh t 



(A.12) 



It would be very interesting to evaluate the functional determinant (|A.ll ) directly and 
check ( |P3 ). 



Appendix B. 



The one-point function ( e iaip ) can be expressed in terms of the appropriately regu- 
larized (see below) Euclidean functional integral 

1(a) = / [VipV$] e~ AM ™ , (B.l) 

JTa 

where Amtm is the Euclidean version of the action (|7|) and the integration is taken over 
the space T a of "twisted" field configurations, such that ip(x) and ip(x) acquire phases 

ip -> e 1 ^ ip, i> $ > ( B - 2 ) 

when continued around the point x = @. In general case it is not known how one can 
evaluate the functional integral (|B.1|) directly. In the case (3 2 — 1/2 the action (^) becomes 
quadratic and the task simplifies drastically. 

The problem is most conveniently treated in conformal polar coordinates (77, 9) 

z = x + iy = e 71 ^ 6 , z = x - iy = e 11 ' 19 . (B.3) 

In this coordinates the action (|7]) with g = takes the form 

A FF = i / dd / ^{^(de-td^L + ^ide + id^R-iMe^ (^ R + ^ L ^ R ) 

JO J -00 

(B.4) 

where ^lr(V:^) an d ^z, r(Vi@) are the components of the Dirac bi-spinors ip and ip 
transformed to the coordinates Q)!. According to (Q) they satisfy the qnasi-periodicity 
conditions 

*l,r(v, + 2n) = _e 2 ^™ ^l,r(v, 0), & L>R (v, + 2n) = _e" 2 ^™ *l,r{v, 0) ■ 

(B.5) 



4 In other words, the fields ip, ip are defined on the universal cover of the punctured Euclidean 
plane R 2 /{0} and satisfy there the quasi-periodicity condition ( |B.2| ). 

5 *Pl{x, y) = e^^ v ^^ e ^ 1,(77, 6), iPr(x, y) = e^^ v+ ^ e ^ r(t], 6), where ipL,R are the components 
of the bi-spinor ip in Cartesian coordinates x,y, and similarly for \td. 
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As usual, the most efficient way to evaluate the functional integral ( |B.1| ) is to use 
the Hamiltonian formalism. The Hamiltonian picture which appears most natural in our 
case is the one where the polar angle 9 is treated as (Euclidean) time, and the Hilbert 
space H is associated with the "equal time" slice 9 = const. We call this approach angular 
quantization®. In this Hamiltonian formalism the fields ty, ^ become operators acting in 
7i. They satisfy canonical equal-time anti-commutation relations 

{t^fa), vj>[(r/)} = 5(n - n') , {y R ( V )^ R (r}')} = 5( V - n') , (B.6) 
while the angular Hamiltonian derived from ( p.4j ) has the form 



K = i 



( dn{-^d v ^ L + ^ R d^ R -Me^ {^ R + ^> L ^ R )] . (B.7) 

J — oo 



The Ty-dependence of the mass term here has the effect of preventing the fermions from 
penetrating too far in the positive rj direction. We call this effect the "mass barrier". 
Correspondingly, the Hamiltonian ( |B.7| ) is diagonalized by the following decompositions 



/OO 1 />00 7 

c uM e~» e , H> R ( V , d)= -= c v u {rj) e~» e , 

roc j ,oo , 

B)= -^=4 <(V) e» 6 , * r (ti, 9) = -= c\ v*M e» e 

J-oo V^TT J-oo V^TT 

in terms of the partial waves 

Ufa) J r(J-tiz) Ui + .(Me^' 1 ^ 



which describe the fermion scattering off the "mass barrier" 
Here 



Angular quantization can be regarded as QFT version of Baxter's Corner Transfer Matrix 



approach |jg|. Angular quantization was previously used in [12| in calculation the sine-Gordon 
form-factors and correlation functions of Jost operators. 
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is the S-matrix associated with this process. In ( |B.8|) c v and c|, are operators satisfying 
the anti-commutation relations 

{c„,c u ,} = {4,cJ,,} = , 

{ci/,4'} = s ( u - v ') ■ 



;b.i2) 



The Hilbert space Ti is the fermionic Fock space over the algebra (P-12|) with the vacuum 
state I v) which satisfies the equations 

c v I v) = for v > , 

4 I v) = for v < . 

The angular Hamiltonian (|B.7| ) can be written as 



(B.13) 



POO 

K = / dv v [clc u + c- v c ] _ v ) 
Jo 



:b.i4) 



As the "Euclidean time" 9 in ( |B.4j ) is compactified the functional integral ( |B.1| ) (with 
g = 0) is given by the trace 



1(a) = Tr n 



-27rK+i^p Q 



(B.15) 



where the the fermion charge 

/OO /"OO 
d77 ( ^ L + tf R V R )= du( 4c„ + c-.cl, ) (B.16) 
-00 Jo 



is introduced in order to impose the twisted boundary conditions Q.2 



The trace (|B.15| ) requires regularization. One can define the system ( |B.4| ) as the 
limiting case e — > of similar system in a semi-infinite box 77 6 [ loge, 00 )0 with the 
boundary conditions 



(B.17) 



Let us denote I e (a) the trace ( |B.15| ) regularized this way. Simple analysis (which takes into 



account the form of the boundary state associated with the boundary conditions (B.17 ) in 
conformal case M = 0) shows that 



(e mv ) = lim e~ 2a J e (o)/J e (0) . 



(B.18) 



7 This corresponds to cutting out a small disc of the size e around the puncture x = in the 



functional integral ( B.l ). The other infinity rj — > +00 does not cause any problem because of the 
"mass barrier". 
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Using the above regularization one can directly evaluate ( [B.15Q , (|B.18|) , with the 



resui 



; e iatp ) = exp 



■ lo § ( ; ,■?. ) d v logS F (v) 



o 



2ni 



(l + e- 2 *")' 



where Sf(v) is the S-matrix ( |B . 1 1| ) . Using the integral representation 



sin(2z/t) 9 . 
— - 2z/e" 2 * 
sinht 



(B.19) 



(B.20) 



one arrives at ( |19|) . 



Eq. ( |i9|) can be checked against exact result of [14J for the spontaneous magnetization 
in the Ising model. Indeed, at the free-fermion point the operator cos(^<£>) can be expressed 
in term of two independent Ising spin fields a^\x) (i = 1, 2) as 



cos[-cp(x 
where a^(x) are normalized as follows 

(T W (X)(7 (1) (0) -> (Xp 1 / 4 



i =2-5 a (1) (x)a (2) (x) 



(B.21) 



Ixl — » . 



According to [|14 



ii i 



:B.22) 



(B.23) 



where A is the Glaisher constant. This constant can be expressed in term of the derivative 
of the Riemann zeta function, 



A = exp( - C'(-l) + ^) = 1.28242712910062... . 



(B.24) 



Evaluating the integral in ([19]) for a = 2 2 and using (|B.24| ) one finds complete agreement 
with (TB723D- 



It is instructive to compare our one-point function ([19]) with known exact results for 
the two-point function fl24] ) in the free-fermion theory ||15|| , ||16|| . For (3 2 = 1/2 and aa' > 
can be written as 



\e iaif (x) e ia,(p (0)) = (e iaifi ){e ia ' lfi } exp( Z aa , (M\x\/2) 



;b.25) 



This result was also obtained by Al. Zamolodchikov [13|. 
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with E aa /(r) denoting the integral 



Y,aa>(r) = - l -J^ pdp j(d pX ) 2 -4sinh 2 X 



,/\2 



2 (a -a' . , 
2 tanhx } ■ 



Here x(p) is the solution of the differential equation 

2 (a — a') 2 

dpX + P~ 1 d P X = 2sinh(2x) H s tanhx ( 1 - tanh 2 x) : 



;B.26) 



(B.27) 



subject to the asymptotic conditions 



X (p) ^V2(a + a') \ log(p) + C aa , 



as p — ► 



-»• ^ sin(v / 27r a) sin(v / 27r a') J K^ {a _ a/) (2p) , as p -> oo 



(B.28) 



and C aa i is a certain constant. Now, one can take the limit \x\ — > in (|B.25|) and use ( |25| ) 
to obtain the relation 



E aa / = Inn ( - 4aa log r + E aa / (r) ) = log ^ ) f^gTT 
where £ a is the one-point function ( |B.19| ) and 



(B.29) 



aaf > 0, |a + a'| < 2 2 . 

We have checked the relation ( |B.29| ) numerically. Note that for pure imaginary a = iu/ y2, 
a' = iu' /y/2 the function x(p) a l so becomes pure imaginary. If in addition u' = u, then 

</>(r) = — 2zx(mr/2) 



satisfies ( |A.5|) , ( |A.7|) in Appendix A, and hence 

If f°° 2 

S(u) = - yE aa ' - A J dp p sinh x 



i=a'=i uj I v2 



(B.30) 



where S(u) is defined by (|A.8|) . The integral in (|B.30| ) can be evaluated by observing that 



4p sinh 2 x = - dp p 2 ( 4 sinh 2 x ~ (d P x) 2 ) 



:b.3^ 



Using (ggg) , ( gig ) one arrives at ( |AJ0|) . 
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Appendix C. 



Expanding ( P3"D in power series in (3 2 one finds 

a 2 = -4 ( 7 + log(ro/2)) + (C(3) - tt 2 ) £ + 0(/3 6 ) , (C.l) 

a 4 = 56C(3)/3 2 + 72C(3)/3 4 + 0(/3 6 ) , (C.2) 

where £(s) is the Riemann zeta function and 7 is the Euler constant. These expansions 
are to be compared with the results of perturbation theory for the action (]5|). In the 
perturbative calculations it is more convenient to use another (equivalent to (^lD) definition 
of the "fully connected" one-point functions a 2n 

a 2 = (^) = lim ( ( <p(z)<p(0) ) I N=e + 4 log e ) , (C.3) 

0~2n = ( (p{xi)(fi{x 2 ) ■ ■ ■ ip(x 2 n) ) c \ Xl=Xa -.- Xan ' ( C - 4 ) 

where (. . .) c in ( |C.4| ) is connected 2n-point correlation function. The Feynman diagrams 
contributing to o 2 and a 4 up to the order (3 A are shown in Fig.l. The calculations are 
straightforward!! and the result coincides with (|C.1|) , ( |C.2|) . 



Expansion of a 2 ( p3|) around /3 2 = 1/2 has the form 

a 2 = -4 (1 + 7 + log(M/2)) + 9 - (7C(3) - 2) + 0(^ 2 ) , (C.5) 



7T 

where (7 = 7r (1/2 — (3 2 )/ (3 2 . The expectation value ( <^> 2 ) can be calculated in perturbation 
theory of the Massive Thirring model ([7]) since (^) allows one to relate the two-point 
function ( (p(x)(p(y) ) in ( |C3| ) to the current-current correlation function 

W(x)j»(0)) = J 0^ (> - U(k 2 ) e ik * , (C.6) 

{(3/2n) 2 ( rtxMO) > = / n ( fc2 ) ^ ■ ( C - 7 ) 

The first two terms of the perturbative expansion in g for the function n(/c 2 ) can be 
obtained by calculating the contributions of the Feynman diagrams in Fig. 2 



M 2 n(/c 5 



1 sinh#-# 



4tt sinh# sinh 2 (#/2) 4tt 2 



(sinh#-0) 2 9 2 



smh z e sinh 2 (#/2) sinh 2 



+ 0(g 2 ) , (C.8) 



One only has to be sure to express the result in terms of the physical mass to. 
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where 6 is related to k 2 as 

k 2 = AM 2 sinh 2 (#/2) . 
Evaluating the Fourier transform ( |C.7| ) and taking the limit ( |C.3| ) one obtains exactly 



Appendix D. 

Exact integral representation for the correlation function 

( vac | of of +n | vac ) 



(d.i; 



( vac | vac ) 

was found recently by Jimbo and Miwa [17]. However in this representation the correla- 
tion function (p.l| ) is expressed in terms of n-fold integral and therefore finding n — > oo 
asymptotic form of ( |L) . 1| ) remains a challenging problem. Here we will explain how the 
asymptotic formula ( |5D| ) , follows from our conjecture (^Dp . 
The XXZ model is a limiting case of XYZ spin chain 

1 oo 

Kxyz = -- ( J * a * <+i + J v ° V s + J z < <+i - J) , (D.2) 

s= — oo 

with J x > J y > | J z |; we introduced here an auxiliary parameter e which will be interpreted 
as a lattice spacing. Following Baxter 0, we use the parameterization 



tt V<9 4 (/5 2 )^i(/3 2 ) # 4 (/3 2 )0 4 (/3 2 

j _ i - p 2 f e^em e 1 (p 2 )e' 1 (o) 



i-p 2 fe' 1 (p 2 ) w 2 ; 



(D.3) 



Jz tt vw 2 ) w 2 )y ' 

J _ l - (3 2 ( e[((3 2 ) 6' 4 ((3 2 ) 



Here 



n V^i(/5 2 ) W 2 ) 

oo 

6»i (u) = 2p* sm(7vv) Yl (1 - p 2n ) (1 - e 2 ™ w p 2n ) (l - e ~ 27riv p 2n ) 

n=l 

oo 

n ( \ TT /'i 2n\ /-i 2ttz u 2n-l\ /i -27ri« 2n— 1\ 

04 ( u ) = \\ (1 -p J(l-e p J (1 — e p J 



n=l 
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and the prime in ( p.3|) means a derivative. Note that for p 
vanishes and the XYZ model reduces to the XXZ chain (13 



the difference J x — J y 



As is known ||18|| , the 



XXZ model (|3l|) describes the critical behavior of XYZ chain, i.e. the correlation length 

i 

R c = eN c — > oo as p — > 0. The scaling limit p — > 0, e ~ p^-e' 2 ) — > of XYZ model is 
described by sine-Gordon QFT, and the following relation between the operators holds 



a x s -> N((3' z ) cos ( ?-(p(r 



ss 



(DA) 



is determined by the constant N, F = \N 2 . 

with the use of two exact 



where N(f3 2 ) is normalization factor. The dependence of n in (|30| ) is of course the simple 
consequence of ( p.4j ) while the factor F in ([ 
It is exact value of this constant which can be deduced from (| 
results in XYZ model: 

a) Exact formula for the energy gap due to Johnson, Krinski and McCoy JHJ, which 
gives in the limit p — > exact relation between p and the physical mass M in the sine- 



Gordon model 18 



», 4 i+e 
M = - p~ 

e 



;d.s) 



where the notation ([TO]) is used. 

b) Exact result of Baxter and Kelland for the expectation value of the spin operator 
a x s [EDI , which in the limit p — > reduces to 



( vac 


a 


X 

s 


vac) 


( vac 


vac) 



(1+0 ps . 



;d.6) 



Comparing (|D.6|) with (|20| ) and using (|D.4|) , (|D.5|) one derives (|32[). 
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